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Abstract 

(-) 

First results on inclusive D and D* ± production in deep inelastic ep scat- 
tering are reported using data collected by the HI experiment at HERA in 1994. 
Differential cross sections are presented for both channels and are found to agree 
well with QCD predictions based on the boson gluon fusion process. A charm 
production cross section for 10 GeV 2 < Q 2 < 100 GeV 2 and 0.01 < y < 0.7 of 
a (ep — > ccX) = (17.4 ± 1.6 ± 1.7 ± 1.4) nb is derived. A first measurement of 
the charm contribution (x, Q 2 ) to the proton structure function for Bjorken 
x between 8 • 10~ 4 and 8 • 10 -3 is presented. In this kinematic range a ratio 
Ff/F 2 = 0.237 ± 0.021+°;°^ is observed. 
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1 Introduction 



The measurement of the charm production cross section in deep inelastic lepton nucleon 
scattering (DIS) is of importance for the understanding of the parton densities in 
the nucleon |]. To order aa s deep inelastic electron nucleon scattering proceeds via 
the generic diagrams depicted in Fig. III. There is some evidence from the EMC 
collaboration || that heavy quark production should be dominated by the boson gluon 
fusion process (Fig. |l]c) 

7 (Z°)S -> cc, (1) 

which has been calculated in next to leading order || . Charm sea quark contributions 
(Fig. 0a,b), where only a single (anti)charm quark is recoiling against the proton, should 
contribute only to a small extent. Charm production due to fragmentation is suppressed 
by heavy quark mass effects and is therefore expected to be small PJ. For the e + e~ 
annihilation process this expectation has been confirmed by recent measurements on 
the Z° resonance 11 |6j . Charmed hadron production via the decay of b flavored hadrons 
is negligible at HERA, because of the small b production cross section. The intrinsic 

a) b) c) 



e e' e e' e e' 




Figure 1: Generic Feynman diagrams for deep inelastic ep scattering up to order aa s : 
(a) ep interaction of the virtual probe (7, Z°) with a valence or a sea quark in the 
proton according to the quark parton model (QPM), (b) corrections to process (a) due 
to gluon radiation off the struck quark before or after interacting with the probe (QCD 
Compton scattering, QCDC), and (c) contribution due to boson gluon fusion (BGF). 

charm model 0, where the proton wave function fluctuates to a \uudcc) state at the 
level of a few permille, will give rise to charm production at large Bjorken x. Currently 
the detection of this process is beyond the scope of HERA experiments. 

A measurement of the charm contribution F^i^x, Q 2 ) to the structure function of 
the proton is supposed to provide information on the nature of the charm production 
process. Present parton density calculations may include heavy flavors in the proton 
via two different routes. One is to include charm in the massless quark evolution || 
starting with a charm density of zero at a scale Qq of the momentum transfer squared 
of the virtual photon, which depends on the mass of the charm quark. This effective 
procedure is presently used in ref. [0, [II], [I^] . The second produces charm exclusively 



via boson gluon fusion taking into account the charm quark mass ||. Although both 
approaches may lead to very similar predictions for F2 C (x,Q 2 ), the resulting charm 
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quark momentum spectra will differ. Therefore the momentum distribution of the 
charmed hadrons allows a sensitive test of the charm production mechanism. 

If the photon gluon fusion process is the dominant source of charmed hadrons in 
deep inelastic ep scattering, the observation of inclusive charmed hadrons production 
represents a sensitive probe of the gluon density xg(x, Q 2 ) in the proton. The measure- 
ment of F2°(x, Q 2 ) can then be used to test different parton density parameterizations. 

(-) 

In this paper a study of inclusive D and D* ± production via the decay channels 

D° -> K-1T+ (2) 

and 

D *+ ^ D \+ ^ K-7V+71+ (3) 

is presented^. After a brief description of the components of the HI detector relevant 
to this analysis (sec. 2) the kinematics of inclusive ep scattering is introduced (sec. 
3). The Monte Carlo simulations used to determine the D meson detection efficiencies 
and to correct the data are briefly described thereafter (sec. 4). The selection of deep 
inelastic charm events and the possible contributions of background sources are then 
discussed (sec. 5). The measured D° and D* + integrated and differential cross sections 
are presented in sec. 6 and compared with QCD based model predictions for heavy 
flavor production in deep inelastic ep scattering. Charm production at HERA is also 
compared with data from vN scattering. The charm contribution, -F| c (x, Q 2 ), to the 
proton structure function is discussed. Finally, the results are summarized in sec. 7. 



2 HI Detector and Data Sample 

In 1994 HERA has been operated with 820 GeV protons colliding with 27.6 GeV elec- 
trons and positrons^. The HI detector [13| is a nearly hermetic multi-purpose apparatus 
built to investigate the inelastic high-energy interactions of electrons and protons at 
HERA. Closest to the interaction point are the central and forward^ tracking systems, 
which are surrounded by a liquid argon calorimeter consisting of an electromagnetic 
and a hadronic section. A super- conducting solenoid surrounding both the tracking 
system and the calorimeter provides a uniform magnetic field of 1.15 T parallel to 
the beam line. The outermost detector component, the instrumented iron, allows the 
measurement of the hadronic energy leaking out of the calorimeter, as well as the 
identification of muon tracks. 

The analysis presented in this paper relies essentially on parts of the central 
tracking system and on the backward electromagnetic calorimeter (BEMC) which are 
described briefly in the following. 

The central jet chamber (CJC) consists of two concentric drift chambers covering 
a polar angle of 15° < G < 165°. It is supplemented by two cylindrical drift chambers 



1 Henceforth, charge conjugate states are always implicitly included. 

2 Subsequently electrons will denote both, electrons and positrons, respectively. 

3 The positive z axis of the HI coordinate system is defined by the proton beam direction. 
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at radii of 18 and 47 cm to determine the z coordinate of the tracks. A cylindrical 
proportional chamber is attached to each of the z drift chambers for triggering. 

In the backward region, the polar angle acceptance of the multi-wire proportional 
chamber (BPC) covers the range 155° < < 174.5°. The reconstructed space point 
together with the z vertex position of the event defines the polar angle of the scattered 
electron with a resolution of 1 mrad. 

The BPC is attached to the BEMC, which covers the polar angular region of 



155° < < 176°. An energy resolution of a(E)/E » 0.1/^jE/GeV © 0.42GeV/£ © 
0.03 is obtained. The energy calibration is known to an accuracy of 1% [[L4j]. A 
scintillator hodoscope (TOF) situated behind the BEMC is used to veto proton induced 
background events based on their early time of arrival compared with particles from 
the interaction at the nominal ep vertex. 

This analysis is restricted to DIS events which have a scattered electron in the 
backward region. These events were triggered by an electromagnetic shower in the 
BEMC with an energy in excess of 4 GeV which was not vetoed by an out of time 
signal in the TOF. The trigger efficiency has been determined from the data using the 
redundancy of the HI trigger system. For energies of the scattered electron candidate 
E' e > 13 GeV, the trigger efficiency is 100%. 

The luminosity is determined from the rate of the Bethe-Heitler reaction ep — > 
epy. The analysis presented in this paper is based on data taken during 1994 run- 
ning periods with electron and positron beams and corresponds to a total integrated 
luminosity of Ci nt = 2.97 pb _1 , with an overall uncertainty of 1.5% |15|. 



3 Kinematics 

At fixed center of mass energy, -^/s, the kinematics of the inclusive scattering process 
ep — > eX is completely determined by two independent Lorentz invariant variables, 
which may be any two of the Bjorken scaling variables x and y, the momentum squared 
Q 2 of the virtual boson and the invariant mass squared W 2 of the hadronic final state. 
In this analysis these variables are determined from the measurement of the energy E' e 
and the polar angle e of the scattered electron according to the expressions (where 
the electron and proton masses are neglected) 



R \ . EL . 2 /0 f 



Q 2 = AE eE ' e cos 2 [^) y = 1 ~Y e sin It 



x = ^- W 2 = Q 2fl X 



(4) 



ys \ x 



where s = 4E e E p , and E e and E p denote the energies of the incoming electron and 
proton, respectively. Here the scattering angle e is defined with respect to the proton 
beam direction. 
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4 Monte Carlo Simulation 



Monte Carlo simulation programs are used in order to correct the data and to estimate 
the systematic uncertainties associated with the measurement. For the determination 
of the acceptance of the detector and the D° and D* + selection efficiencies, heavy flavor 
(charm and bottom) DIS events are generated using the AROMA 2.1 |Tj| program. 
This simulates neutral current heavy quark production via the boson gluon fusion 
process, which is implemented in leading order QCD including heavy quark mass effects. 
A charm quark mass of m c = 1.5 GeV is used. Higher order QCD radiation includes 
initial and final state parton showers (PS) [0 using the LEPTO 6.1 code fll8 |. JETSET 



7.4 |]T9| is used to perform the hadronization according to the Lund string model with 
the symmetric Lund fragmentation function | 2"0fl . For the gluon density in the proton 
the MRSH parameterization [] f is used. This is compatible with measurements of the 
proton structure function at HERA [21]. In order to test the sensitivity of the results to 
the charm quark mass heavy flavor DIS events are also generated with m c = 1.3 GeV 
and m c = 1.7 GeV. 

For the study of any charm contribution to the proton sea, heavy flavor events are 
also generated using the LEPTO 6.1 program In this generator heavy quark mass 
effects are not included. To avoid the divergences in the QCD matrix elements, a cut 
in the smallest invariant mass rriij of any two partons of an event is introduced. The 
difference in the cross section, obtained with this cut, to the total DIS cross section 
is then attributed to sea quarks in the proton, thus also generating QPM type charm 
events, where a single charm quark recoils against the proton remnant. 



5 Event Selection 

5.1 Selection of Deep Inelastic Scattering Events 

For the selection of deep inelastic ep interactions the electron candidate is identified 
as the particle giving in the BEMC the most energetic cluster, which has also to fulfill 
the following conditions : (a) its center of gravity has a radial distance of less than 4 
cm from a reconstructed BPC point, and (b) its lateral shower size has to be smaller 
than 4 cm. In order to calculate the kinematic variables precisely the position of the 
event vertex is needed. It is defined by at least one well measured track crossing the 
beam axis. 

The basic kinematic constraints are that the electron polar angle is in the range 
155° < © e < 173° and that the electron energy E' e > 13 GeV. The latter requirement 
ensures that the trigger is fully efficient and the photoproduction background is small 

The analysis is thus restricted to 10 GeV 2 < Q 2 < 100 GeV 2 and y < 0.53. 

The dominant source of non-ep background in the event sample is due to interac- 
tions of the proton beam with the residual gas and the beam collimators before the HI 
detector. The level of this background is studied with events originating from proton 
bunches having no colliding electron bunch partner. It is found to be less than 0.25% 
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Figure 2: Distribution of (a) the polar angle and (b) the energy of the scattered electron 
for charm events as observed in the D* + analysis. The background subtracted data is 
compared with the expectation of the AROMA (histogram) Monte Carlo simulation, 
normalized to the number of events in the data. Only statistical errors are shown. 



of the total number of the selected events. 

Events originating from photoproduction, in which the scattered electron escapes 
detection but hadronic activity in the BEMC fakes an electron, are the only significant 
ep related background. On the basis of Monte Carlo simulations, using the PYTHIA 
generator [Tj|, this background is found to be less than 0.4%. 



Figure [| shows the distribution of the angle and the energy of the scattered elec- 
tron for charm events as selected by the D* + analysis in comparison to the expectation 
of the AROMA Monte Carlo simulation. Good agreement is observed between the 
background subtracted data and the simulation. 



5.2 Selection of Charm Events 

Reconstructed D°(1864) and D* + (2010) mesons are used to tag heavy quark produc- 
tion. The analysis is based on charged particles reconstructed in the CJC. Apart from 
the electron candidate at least 2 or 3 tracks, depending on the decay mode, have to 
be fitted to a common event vertex. No particle identification is applied so each track 
may be assumed to be a pion or a kaon. A minimum track length of 10 cm in the 
plane perpendicular to the beam direction is required for those tracks considered to 
form the D° candidate. In both analyses the K~ir + mass combination has to lie within 
the interval \r] K7T \ < 1.5, where r\ = — In tan 0/2 denotes the pseudorapidity. 

The D° is identified via its decay mode 

D° -> K-7i + (5) 
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and the D* + through the decay chain 

D*+ -> D\f low -> tt+tt+u, . (6) 

For the latter use is made of the tight kinematic constraint for the decay of D* + — > 
D n^ low [p2| . A better resolution is expected for the mass difference 



Am = m(D°7r+ low ) - m(D°) (7) 

than for the _D* + mass itself, the width of which is dominated by the momentum 
resolution of the detector. No cut on the track length is applied in selecting the slow 
pion. 

In addition to the reconstruction of the specific final states the general strategy 
in searching for charm induced events is based on the hard fragmentation of charmed 



hadrons p3 |. The analysis is performed in the hadronic center of mass (7*2?) system. 



In analogy to the analyses performed by e + e experiments the quantity 



x D 





- 2 


Pdo 






W 



is defined for the D° in the 7*p frame for both decay channels^. The mean (xd) of 
particle combinations originating from D° mesons is found to be significantly larger 
than for the combinatorial background. It therefore may be expected that the D° 
decay products have large momenta in the 7*p system, especially in view of the low 
multiplicity of the decay mode used in this analysis. In order to take advantage of 
this fact a ranking of the charged particles is introduced, such that the particle of a 
given charge having the largest momentum in the 7*p system gets the rank R = 1, the 
particle with the second largest momentum gets the rank R = 2, etc. Cut are applied 
in both, the normalized momentum xd and the rank of the particles in the subsequent 
analyses in order to suppress the combinatorial background. 



5.2.1 Selection of D* + -> D°7r+ ow -> K~tx + ti+ ow Decays 

Only particles having a transverse momentum of p t > 0.25 GeV in the laboratory frame 
are considered for the reconstruction of D° candidates. The mass combinations 
fulfilling the requirements 

(R K = 1 A R n < 3) V (R w = 1 A R K < 2) (9) 

are calculated for each event. The asymmetric cut in the order of the kaon (Rk) and the 
pion (R n ) accounts for the effect of their different masses. The fractional momentum of 
the mass combinations has to satisfy xd > 0.25. At small Bjorken x large xd favors D° 
candidates going backward in the laboratory frame. Also those irk, mass combinations 
are accepted, which satisfy the condition p t (K~7i + ) > 3 GeV, irrespectively of Rk, R-k 
and xd, to retain the D* + mesons going in the forward direction at small Bjorken 
x. The rriKn combinations within ±90MeV of the nominal D° mass of 1.865 GeV are 
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Figure 3: Distribution of (a) the mass difference Am and (b) the K~7t + mass for DIS 
events selected as described in the text. The data points in (a) are obtained from the 
K~7i + mass combinations fulfilling \rriKn — m D°\ < 90 MeV. The shaded histogram 
shows the background expectation, which is obtained from the high mass sideband 
2.05 GeV < m Kn < 2.5 GeV normalized to the region 0.160 GeV < Am < 0.180 GeV. 
The data points in (b) are obtained from the candidates in an interval of 2.2 MeV 
around the nominal Am value expected for the decay D* + — > D°tt + . The shaded 
histogram shows the background expectation from the region 0.170 GeV < Am < 
0.180 GeV normalized by the two body phase space factors according to the Am 
intervals. The solid lines represent the result of the fits as described in the text. 

combined with each additional charged track with p t > 0.12 GeV and a charge opposite 
to that of the kaon candidate. 

Figure |^a shows the distribution of the mass difference for accepted K~7t + tt + 
combinations, fulfilling the requirement [mjcn — m D o\ < 90 MeV. Clear evidence for 
D* + production is seen in the signal region, defined by a ±2.2 MeV window around 
the expected mass difference of mr>*+ — m^o = 145.4 MeV. The fit to the data includes 
the D* + signal, partially reconstructed D° mesons and the decay modes D° — > 7r + 7r~ 
and D° — > K + K~~ , and the two body phase space background. An inclusive D* + yield 
of 103 ± 13 is observed. The peak position at Am = 145.5 ± 0.15 MeV as well as 
the width of the Gaussian of a = (1.04 ± 0.18) MeV agree well with the Monte Carlo 
expectations. No enhancement is observed in the background sample obtained from 
mass combinations with 2.05 GeV < mKn < 2.5 GeV as shown by the shaded histogram 
in Fig. |Q. 

Figure|3|b shows the mass distribution of the K~tt + combinations for events in the 
Am signal region. The shaded histogram shows the background expectation from the 
region 0.160 GeV < Am < 0.180 GeV scaled by the two body phase space factors. It 

4 Quantities defined in the 7*p system are marked by *. 

5 Since the selection is mainly based on the leading particles in the 7*p frame, the background 
behavior cannot be studied by using the like sign K ~n + mass combinations. 
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Figure 4: K~ix + mass distribution observed in the inclusive D° analysis. The solid line 
is a fit to the data as indicated in the text. The dashed line shows the contribution from 
the wrong K — 7r assignment to the tracks. The shaded area indicates contributions to 
the fit due to combinatorial background and reflections while the dashed line describes 
the mass distribution for the wrong mass assignment to the tracks. 

describes well the K~n + mass distribution for > 2.05 GeV. Below 2.05 GeV large 
differences are observed in the rriKw distribution from the D* + signal region as compared 
to the background expectation. These differences are due to contributions where both 
particles originate from the decay D* + — ► D°7if low . Clear signals are observed for 
the decays D° — > K~rr + and D° — ► K~tt + 7t° where the ir° is not detected. The fit 
to the data takes into account these two contributions over an exponentially falling 
background. The inclusive D° yield of 84 ± 16 for the decay D° — > K~rr + agrees with 
the D* + rate obtained from the fit to the Am distribution. 

5.2.2 Selection of D° -> K~n + Decays 

For each event the K~ir + mass combination is calculated only for those particles with 
p t > 1 GeV and fulfilling 

R K = 1 A R n < 2 . (10) 

For the K~tt + mass combination a p t > 2 GeV and xd > 0.3 is required. Finally a 
cut in the helicity angle 0^- of the kaon in the D° rest frame with respect to the D° 
direction in the laboratory system of | cos6^| < 0.5 is imposed. 

Figure [| shows the mass distribution of the K~n + combinations accepted after 
these cuts. Clear evidence for D° production is observed. The fit to the data contains 
contributions for the D° with the correct K, it assignment to the tracks, for the D° 
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with the wrong K, ir assignment to the tracks, an exponential distribution for the 
combinatorial background, and contributions due to reflections from other D° decay 
modes and from D + decay at small m^. An inclusive D° signal of 144 ±19 is observed 
within the cuts. In total 20 ± 5 events are found to be selected by both the D° and 
D* + analyses. 



6 Results 

The results on the differential and the integrated cross sections for the production of 
charmed hadrons in ep collisions are presented for the range 10 GeV 2 < Q 2 < 100 GeV 2 
and 0.01 < y < 0.7 with a small, 7.8% contribution of the extrapolation from the 
experimentally accessible range of y < 0.53 to 0.7. The cross sections are calculated 
from the observed numbers iV & s of D° and D* + candidates, according to 

a(ep^DX)= % , (11) 

>->int ' a " ^tot 

where a(ep — > DX) is used as an abbreviation for both the differential and integral 
cross sections, respectively. D stands for D° and D* + and their charge conjugates. 
Here Ci n t and B refer to the integrated luminosity and the branching ratios B = 
B(D° -> K-Ti + ) = 0.0401 ± 0.0014 and B = B(D* + -> D°ir+) ■ B(D° -> K'^) = 
0.0273 ± 0.0011 p3| . The quantity e to t is the total efficiency. 



The number of observed events, N f, s , is obtained from the fits to the spectra. In 
the case of differential distributions the mass distribution for each bin is fit separately 
by fixing the position and the width of the signals to those values obtained from the 
total data sample. The effect of correlations between the detector resolution and the 
measured quantity is included in the experimental systematic errors. The AROMA 
Monte Carlo program is used to determine the efficiency arising from the selection 
procedure and the geometrical acceptance of the detector. A total efficiency, including 
the acceptance of the apparatus, of (16.3^2^)% and (5.9^ ' 6 )% is obtained for the 
inclusive D* + and the inclusive D° analyses, respectively. 

The inclusive D° and D* + cross sections are then converted into a charm produc- 
tion cross section in deep inelastic ep scattering using the equation 

<j(ep — * eccX) = - ■ — ; — ; , 12 

v y 2 P(c —>£>)• (1 + £) V ; 

where P(c — > D) denotes the charm quark fragmentation probability into a specific D 
meson and £ the correction to apply for charm production via fragmentation or via B 
meson decays, respectively. 

The fragmentation probabilities P(c — > D°) and P(c — > D* + ) have to be deter- 



mined from other experiments. The ARGUS [|24| and CLEO |25| results from below the 
T(4S*) together with the HRS measurement assuming P(b — > D°) = P(c — > D°) 



are averaged to get P(c -> D°) ■ BR(D° -> K~7r + ) = 0.0205 ± 0.0011. For L> H 
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a(ep -> eM) [nb] 


cr(ep — > eccX) [nb] 


D° — ► 

D* + - 


K-TT+ 

* (K-lT + )lT + 


20.4 ± 2.7 -il 
7.8 ± 1.0+J;J ±0.6 


i9.5 ± 2.6 + 2 : ™ 

15.1 ± 1-8+2^ ± 1.2 


Averaj 


Z,C 




17.4 ±1.6 ±1.7 ±1.4 



Table 1: Inclusive D meson and charm production cross sections in deep inelastic ep 
scattering for the kinematic range 10 GeV 2 < Q 2 < 100 GeV 2 and 0.01 < y < 0.7. 
The errors refer to the statistical error, the experimental systematic, and the model 
dependent uncertainties. 

mesons the average of all e + e~ data below the Z° leads to P(c — ► D* + ) ■ BR(D* + — > 
D°tt + ) ■ BR(D° -> R-7T+) = 0.0069 ± 0.0005 according to the procedure of ref. HQ. 

The main contribution to the correction £ arises from gluon splitting. Its esti- 
mation is based on the OPAL measurement of the mean multiplicity per event for cc 
creation due to gluon splitting with a value (n g ^ C c) = 0.0238 ± 0.0048 0. This result 
has been extrapolated to HERA energies yielding £ s ^ C c = 0.02 ± 0.02 for xp > 0.25. 



Experiment 


ref 


Process 


CT (D*±X) 
a(D°X) 


HI 






ep 


0.38 ± 0.07 










±0.06 


CLEO 




25 


e + e~ 


0.48 ±0.06 


ARGUS 




24 


e + e~ 


0.47 ±0.06 


HRS 




26 


e + e~ 


0.47 ±0.06 


DELPHI 




27 


e + e~ 


0.43 ±0.05 


ALEPH 




28 


e+e" 


0.36 ±0.04 


OPAL 


ii 


| e+e" 


0.38 ±0.03 


NA32 




m 


7T Si 


0.54 ±0.05 


E769 




31 


n Be,Cu,Al,W 


0.39 ±0.05 








p Be,Cu,Al,W 


0.32 ±0.13 








K Be,Cu,Al,W 


0.24 ±0.08 



Table 2: Experimental Ratio of the inclusive D*^ cross section over the inclusive D 
cross section compared to the results from e + e~ and hadroproduction experiments. 
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Model 


m c [GeV] 


Predictions [nb] 


This Experiment [nb] 


GRV 


1.5 


11.4 


17.2 ± 2.3 


MRSH 


1.3 


11.3 


17.2 ±2.3 


MRSH 


1.5 


9.7 


17.4 ±2.3 


MRSH 


1.7 


8.1 


16.7 ± 2.2 


MRSA' 


1.5 


9.7 


17.4 ±2.3 


MRSDO' 


1.5 


8.7 


18.3 ±2.4 


HI F 2 fit 


1.5 


13.6 ± 1.0 


17.1 ±2.3 



Table 3: Comparison of predicted charm production cross sections with the present 
measurement for 10 GeV 2 < Q 2 < 100 GeV 2 and 0.01 < y < 0.7. The theoretical 
predictions are based on NLO calculations for different parton densities and different 
values of m c . Also shown is the prediction based on the determination of the gluon 
density by the HI NLO fit to the total F 2 data ||15|| . The experimental cross section is 



the average from the D* + and D analysis. Total experimental errors are given. 



6.1 Integrated Cross Sections 

In Tab. [I] the inclusive D meson and charm production cross sections at the Born level 
are summarized for the kinematic range 10 GeV 2 < Q 2 < 100 GeV 2 and 0.01 < y < 0.7. 
Corrections for QED radiations are applied using the HECTOR program [jSHfl . They 
amount to 5 ± 3 % for the kinematic region explored in this analysis. The errors 
in Tab. [1] reflect the statistical uncertainties, the experimental systematic and the 
model dependent uncertainties due to the extrapolation to the full phase space. The 
experimental systematic errors and the model dependent uncertainties are discussed 
separately in the subsequent sections. The results of both measurements are in fair 
agreement. Table p] also includes the charm production cross section obtained by 
combining the results of the two analyses. 

The ratio on the inclusive D* + over the inclusive D° cross section is compared with 
results from e + e~ annihilation and hadroproduction experiments in Tab. |2|. Agreement 
is observed with the other experiment s[]. The quoted errors for the present measure- 
ment refer to the statistical and the experimental systematic error, and take into ac- 
count the systematic uncertainties, common to both analyses. 

In Tab. |3] the measured charm production cross section is compared to the NLO 
predictions [[3] for different parameterizations of the gluon density in the proton and 
different values of the charm quark mass m c . The results from the inclusive D° and 
D* + analyses are combined for the determination of the charm production cross section. 
The quoted error refers to the total experimental error, which is obtained by adding 

6 The value quoted here differs slightly from that given in Q because it also includes the measure- 
ment of ARGUS H . 

7 The ratio from the ARGUS and HRS experiment is calculated using a value of i?i?(_D* ± — * 
D°7r±) = 0.681 H p. 1171] instead of 0.55 as quoted in their paper [|| |§. 
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Experimental 



Systematic Errors 


D° 


D*+ 


rpflpr't ion d 






hrirkpTonnd shanpi^ 






mass resolution 


±0.10 


±0.11 


wrong if _ 7r + assignment 


±0.03 


- 


tracker efficiency 


+0.06 
-0.02 


+0.09 
-0.03 


luminosity 


0.015 


0.015 


7P contribution 


<0.004 


<0.004 


radiative corrections 


±0.03 


±0.03 


branching ratios 


±0.03 


±0.04 


total ep — > eDX 


+0.13 
-0.12 


+0.15 
-0.13 



Table 4: Summary of the relative experimental systematic uncertainties of the inte- 
grated inclusive cross section. 



the statistical and the experimental systematic errors quadratically. The measured 
cross section is fairly insensitive to the choice of the model assumptions. For the 
theoretical predictions, however, the integrated charm production cross section shows 
strong variations on the model assumptions. All theoretical predictions based on the 
NLO calculations |§ are below the measured cross section. 

The table also includes the comparison of the data with the expectation of the 
NLO QCD analysis of the structure function data performed by the HI Collabora- 
tion |15[. In the evolution only three light quark flavors are taken into account. Heavy 
quark contributions are dynamically generated using the BGF prescription given in 
ref. P3|, extended to NLO according to ref. H. The scale of the BGF process 



has been taken as yQ 2 + 4m 2 with a charm quark mass of m c = 1.5 GeV. The gluon, 
sea and valence quark distributions were parameterized at Ql = 5 GeV 2 . The QCD 
parameter Aq C d is kept fixed to the value of 263 MeV, as determined in ref. PS] . The 



parton densities are derived from a fit of the evolution equations to the data [15, 36, 37 



The error on the predicted charm production cross section was obtained in propagat- 
ing the statistical and the uncorrelated experimental systematic errors on the total F 2 
data through this fitting procedure [pSfl . This determination of the gluon density in 
the proton yields a predicted charm production cross section which comes closer to the 
result of the current analysis. 



6.1.1 Experimental Systematic Uncertainties 

The experimental systematic uncertainties are summarized in Tab. f|. By adding the 
different contributions in quadrature total relative errors of and are found 
for the inclusive D° and D* + cross sections, respectively. 
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Model 



Uncertainties 


D° 


D* + 


charm fragmentation 






function 


±0.04 


±0.04 


GRV-HO Gluon 


-0.01 


-0.01 


MRSDO' Gluon 


+0.05 


±0.05 


Gluon from HI F 2 fit 


-0.01 


-0.02 


m c = 1.3 GeV 


±0.04 


±0.03 


m c = 1.7 GeV 


±0.04 


-0.06 


total ep — > eDX 


+0.08 
-0.06 


±0.07 


fragmentation & 






gluon splitting 


±0.02 


±0.02 


6 contribution 


<0.02 


<0.02 


total ep — > eccX 


+0.08 
-0.06 


±0.08 



Table 5: Summary of the relative model dependent uncertainties of the inclusive D 
meson and charm production cross sections. 



In the case of the inclusive D° analysis the dominant experimental systematic 
uncertainty arises from the description of the combinatorial background and the shape 
of the reflections below the D° mass. A simple phenomenological ansatz has been 
used to parameterize these contributions in the fit shown in Fig. [|. The uncertainty 
introduced by this procedure is estimated by varying the txik-k range in which the data 
is fitted. The form of the wrong K — n assignment and its rate with respect to the 
correct K — tc assignment is studied by Monte Carlo simulation. The given error is 
determined from this simulation. 

For the inclusive D* analysis, the uncertainties coming from the shape of the 
combinatorial background are determined from the Am distribution of the high m& 
sideband, which matches well the two body phase space expectation. The given uncer- 
tainty accounts for the fact that the shap of the combinatorial background may differ 
from the two body phase space prediction. The contribution due to reflections in the 
D° mass window is studied with Monte Carlo methods. 

For the differential cross sections in the subsequent sections an error of 10% due 
to the procedure used to determine the inclusive D meson yields per bin is added 
quadratically to the systematic uncertainties. 

The uncertainty of the single track reconstruction efficiency are investigated as 
described in ref. |3TI . 
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6.1.2 Model Dependent Uncertainties 



Various sources of model dependent uncertainties of the cross sections have been inves- 
tigated and are summarized in Tab. |5|. They are based on estimated errors on the total 
efficiency due to the uncertainties in the charm fragmentation function and its QCD 
evolution, in the gluon density of the proton, and in the effect of the charm quark mass 
on the inclusive D cross section. In the case of the charm production cross section, 
they also include an error for the estimate of the correction factor £. The total model 
dependent error is obtained by adding the different contributions quadratically. 

The influence of the charm fragmentation function on the cross section is checked 
by using different parameterizations |D|. A change in the total efficiency of less 
than 1% is found. The effect of the QCD evolution is studied by shifting the mean (xd) 
by 0.1 which corresponds roughly to the difference in the gluon radiation by changing 
W by more than one order of magnitude. This results in an uncertainty of 4%. 

The sensitivity of the measured cross section on the parton density in the pro- 
ton has been studied by using different parameterization. In addition to the MRSH, 
MRS A', MRSDO' and GRV-HO parameterizations also the gluon density resulting from 
the NLO QCD fit to the total F2 measurement from HI ]T3| has been used for efficiency 
calculations. Only small variations in the total efficiency are observed. 

In order to investigate the influence of the charm quark mass the correction 
factors have been calculated by using the AROMA Monte Carlo program with m c = 
1.3 GeV, 1.5 GeV, and 1.7 GeV. Again only small variations in the efficiency are 
observed. 



6.2 Differential Cross Section 

Figure [5| shows the normalized differential cross section 1/cr da / dp t (D°) for inclusive 
D° and for D° via D* + production for 10 GeV 2 < Q 2 < 100 GeV 2 and 0.01 < y < 
0.7. The normalization is done with respect to the integrated cross section. Only 
statistical errors are shown. No source of systematic uncertainties is found which 
would introduce additional uncorrelated systematic errors. Good agreement between 
the normalized pt(D°) spectra of the two analyses is observed^. The data are also 
compared with the expectation from the boson gluon fusion process for charm quark 
masses of m c = 1.3 GeV and m c = 1.7 GeV. The measured differential distributions 
are in good agreement with the expectation. The differences due to the choice of the 
charm quark mass are much smaller than the present accuracy of the data. The data 
have also been compared to calculations using the GRV-HO |§ and MRSDO' ]10|] gluon 



densities, which lead to normalized pt(D°) distributions nearly identical with the one 
obtained for the MRSH parameterization. 



8 Although the D° produced via D* + decay does not originate directly from the charmed quark, 
no differences between the pt(D°) spectra for the two analyses are expected due the smallness of the 
kinetic energy available in the D* + — > D°tt + decay. 
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Figure 5: Normalized distributions of 1/a da / dp t (D°) of the D° mesons from the 
D* + (full points) and the D° analysis (open circles) in comparison to the AROMA 
expectation for m c = 1.3 GeV (full line) and m c = 1.7 GeV (dashed line) using the 
MRSH parameterization of the parton densities in the proton. 



6.3 Mechanisms of Charm Production in ep Scattering 

The charm production mechanisms in ep-scattering can be investigated by comparing 
fully corrected data with expectations for charm production off a charm quark in the 
proton sea and via boson gluon fusion (BGF). It is assumed that the fragmentation 
function D c D (z) with z = Ed/E c , the ratio of the charmed meson energy to the 
charm quark energy E c , is a universal function independent of the charm creation 



mechanism. This has been shown in ref. [^T| to hold for both electron-positron an- 
nihilation and neutrino-nucleon scattering [42, (E| provided the data is corrected for 
QCD radiation ||. Figure^ shows the corrected distributions 1/a da/dx^, which are 
a convolution of the charm production spectrum with the fragmentation function in 
comparison with: 



1. Monte Carlo simulated data using the AROMA generator in which charmed 
quarks are produced by BGF and fragmented according to JETSET. 

2. Measurements of V iV-scattering by the CDHS jEJ and E531 JEJ experiments, 
where charm production is expected to proceed mainly via W /± -scattering off 
a strange sea quark in the proton^ (in analogy to Fig. la). It can be shown 
that the variable z, measured in the proton rest frame, approximately transforms 

9 In case of the vN scattering, the data also include sizeable contributions from valence quarks. 
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Figure 6: Normalized xd distribution in deep inelastic ep scattering at (W) ~ 125 GeV 
for \t]u\ < 1.5. The open points represent the D° data, the closed points the D* 
data. The shaded histogram shows the expectation of the boson gluon fusion process 
according to the AROMA Monte Carlo simulation. The dashed histogram shows the 
expectation from a charm sea contribution, calculated by selecting QPM events with 
the LEPTO/MEPS Monte Carlo program. The full line gives the result of the QCD 
evolution of the v N data. 

into xp in the W*p system. The results are evolved from (W) = 8 GeV to 
(W) = 125 GeV ||. 

3. Monte Carlo simulated data using the LEPTO 6.1 generator, from which only 
charm sea quark events are selected. The fragmentation is modelled using JET- 
SET. 

The data are restricted to the region \t]d\ < 1-5 in order to avoid large xd 
dependent correction factors due to the detector acceptance. The total error is obtained 
by adding in quadrature the statistical error for the real and simulated data and the 
sum of the model dependent systematic uncertainties. The measurement errors are 
dominated by statistical effects. No source of systematic uncertainties is found which 
would introduce significant uncorrelated errors. Good agreement is observed between 
the xd distribution obtained from the two different analyses. 

The differences between the expectation from the BGF model, where two charm quarks 
are recoiling against the proton remnant in the hadronic center of mass system, on the 
one hand and the predictions from the v N data and the quark parton model, where 
only one charm quark emerges opposite to the proton remnant, on the other hand are 
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Mode (Q 2 ) (x) Ff Ff/F 2 

[GeV] 



D* 


12 


.0008 


0.211 ±0.049 


+0.045 
-0.040 


0.198 ±0.045 


+0.039 
-0.033 


D° 


12 


.0020 


0.263 ±0.036 


+0.043 
-0.041 


0.304 ±0.043 


+0.045 
-0.042 


D* 


12 


.0032 


0.190 ±0.054 


+0.052 
-0.049 


0.254 ±0.072 


+0.063 
-0.058 



D* 


25 


.0008 


0.324 ±0.099 


+0.065 
-0.058 


0.244 ± 0.075 


+0.044 
-0.038 


D° 


25 


.0020 


0.253 ±0.069 


+0.041 
-0.040 


0.248 ± 0.086 


+0.035 
-0.033 


D* 


25 


.0032 


0.222 ±0.066 


+0.044 
-0.039 


0.255 ±0.076 


+0.045 
-0.038 



D* 


45 


.0020 


0.156 ±0.070 


+0.031 
-0.028 


0.127 ±0.059 


+0.023 
-0.020 


D° 


45 


.0032 


0.275 ± 0.074 


+0.045 
-0.042 


0.249 ±0.071 


+0.035 
-0.033 


D* 


45 


.0080 


0.200 ±0.064 


+0.040 
-0.035 


0.269 ±0.088 


+0.047 
-0.039 



Table 6: The value of F£ c and of the ratio F| c /F 2 . The errors refer to the statistical 
and the experimental systematic errors. 

evident. The BGF model agrees very well with the data which justifies the usage of 
this model throughout the paper. 

To determine the fraction of the charm production cross section which may be 
attributed to a photon scattering off charm sea quarks in the proton, a superposition 
of the boson gluon fusion and the sea quark predictions, i.e. a function of the shape 

1/ada/dxo = (1 — e)(l/ ada / <Ixd)bgf + e (l/ 'vda j l dxri) sea (13) 

is fitted to the xd distribution of Fig. ^. Depending on the form of a charm sea quark 
contribution, i.e. using the results of the LEPTO Monte Carlo generator or the V N 
data, and depending on the charm quark mass m c , the fit yields values for e between 
—0.062 ± 0.035 and —0.041 ± 0.031. According to the Bayesian approach ||23] these 



values are converted to an upper limit on a charm sea quark contribution to charm 
production in deep inelastic ep scattering in the kinematic range at HERA of 

e < 0.05 (14) 

at the 95% confidence level. It may thus be assumed that boson gluon fusion is the 
dominant charm production process in DIS at HERA. 

6.4 Charm Contribution to the Proton Structure Function 

The charm contribution -F| c (x, Q 2 ) to the structure function is obtained by using the 
expression for the one photon exchange cross section for charm production 

(l±(l-2/) 2 ) F 2 cc (x,Q 2 ), (15) 



dxdQ 2 Q 4 x 
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Systematic Errors 
on F£ c 


D 


D*+ 


fit 


0.14 


0.15 


track efficiency 


+0.06 
—0.02 


+0.09 
—0.03 


bin miPT^.tinns 


0.02 


0.05 


1 LLlillllUol by 




VJ.VJIO 


7P contribution 


<0.004 


<0.004 


radiativp porrpptiorm 


0.03 


0.03 


ii agmen i a i ion 






gluon splitting 


0.02 


0.02 


6 contribution 


<0.02 


<0.02 


charm fragmentation 






function 


0.04 


0.04 


P(c -> D) 


0.05 


0.07 


total 


+0.17 
-0.16 


+0.20 
-0.19 



Table 7: Summary of the average relative systematic uncertainties on the measurement 
of Ff. 

with the simplification that the Callan-Gross relation holds, i.e. R = F 2 /2xFx — 1 = 0. 
According to QCD calculations a maximum value of R ~ 0.38 is expected in the 
kinematic range explored in the present analysis. An increase in F^^x, Q 2 ) of at most 
2% is obtained at small x for R — 1. 

The charm contribution to the proton structure function is obtained from the 
numbers of reconstructed D* + and D° mesons, which are converted to bin averaged 
cross sections according to Eqns. ( |Il~| , |T2"D and using the Monte Carlo efficiency calcu- 
lation. The coarse binning in x and Q 2 is dictated by the small statistics available. 
The bin averaged cross section is corrected for QED radiative effects using the program 
HECTOR [0. 

In Tab. || the results of the F| c measurements are summarized. The errors refer to 
the statistical and to the experimental systematic uncertainties. The systematic uncer- 
tainties are summarized in Tab. 0. They include an error due to the angular and energy 
resolution of the electron measurement, which is dominated by the uncertainty of the 
BEMC energy calibration. A shift in the energy scale of 1% introduces uncertainties 
of 6% (9%) in Ff for the low (high) x bin in the D* + analysis at (Q 2 ) = 12 GeV 2 due 
to bin migrations in x. For all other bins this effect is below 5%. The experimental 
systematic errors also include the uncertainties in the determination of P(c — > D* + ) 
and P(c^ D°). 

The measurement of _F 2 CC is found to be fairly insensitive to the actual choice of 
the parameterization of the gluon density in the proton. Using GRV-HO or MRSD0' 
parameterization of the parton densities instead of MRSH, changes the result by less 
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Figure 7: The charm contribution F 2 CC to the proton structure function as derived 
from the inclusive D* + (full dots) and D° analysis (open circles) in comparison with 
the NLO calculations based on GRV-HO (full line), MRSH (dashed line), and MRSDO' 
(dash-dotted line) parton distributions using a charm quark mass of m c = 1.5 GeV 
for (Q 2 ) = 12, 25 and 45 GeV 2 . The inner (outer) error bars refer to the statistical 
(total) errors. The shaded band represents the prediction from the HI NLO fit to the 
F2 measurements. The EMC data are also shown (open boxes). 



than 2% and 5%, respectively. 

Table M also includes the ratio FSfjFi where the Fi{x, Q 2 ) measurement is taken 



from ref. H15] . Within the accuracy of the charm data this ratio seems to be constant 
in the kinematic range explored in the current analysis. In particular the steep rise of 
F^jF% with Q 2 for fixed x, as found by the EMC experiment 0] at considerably larger 
x, is not observed at HERA. The average yields 

+0.043 
-0.039 



f/F 2 ) = 0.237 ± 0.021 " , (16) 



which is an increase of about one order of magnitude of the overall charm contribution 
compared to the EMC result. This is consistent with the measured rise of the gluon 
distribution towards low x and the dominance of the BGF process observed here. The 



present result on F 2 CC / ' F2 is also found to be in range predicted by ref. p4 

The F% c measurements are displayed in Fig. |7| together with the result of the 
EMC collaboration ||. The measurement at HERA extends the range of the _F 2 C 
measurement by two orders of magnitude towards smaller x values. The comparison 
of the HI and EMC measurements reveals a steep rise of _F 2 CC with decreasing x. The 
data are compared with NLO calculations |J using the GRV-HO, the MRSH, and 
the MRSDO' parameterizations of the gluon density in the proton for a charm quark 
mass m c =1.5 GeV. So far no experimental information is available on the gluon 
density in the proton observed in charm production at small x. Therefore, also the 
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MRSDO' parameterization is considered here, although it has already been excluded 
by inclusive measurements at HERA f45fl . This parameterization also fails to describe 
F 2 CC . The data lie systematically above all predictions which is to be expected from the 
measurement of the charm production cross section. 

The data are also compared to the prediction from the HI QCD fit to the F 2 
measurements using a charm quark mass of m c = 1.5 GeV. The error band shown for 
this fit includes the propagation of the statistical and the uncorrelated systematic errors 
on the total F 2 data through the fitting procedure. This prediction is systematically 
above all other calculations, independently of x and Q 2 but agrees better with the 
F 2 C measurement. The charm data indicate a small excess of the measured F 2 C with 
respect to all predictions in the lowest Q 2 bin. 

The dominant uncertainty in the QCD calculations arises from the uncertainty in 
the charm quark mass. This affects mainly the lowest Q 2 bin, for which a variation of 
m c by 200 MeV will change the prediction for F 2 C by 15%. However, the measurement 
of F 2 C shows the same m c dependence. Therefore the comparison of the data with QCD 
calculations based on different gluon densities in the proton is nearly independent on 
the assumption made for the charm quark mass. 



7 Conclusions 

Results on inclusive D° and D* + meson production in neutral current deep inelastic 
ep scattering at HERA have been presented. It has been shown, that the production 
dynamics of charmed mesons in the current and central fragmentation region may be 
described by the boson gluon fusion process. 

The observed inclusive cross section ratio a[ep — > D* + X) / <j(ep — > D°X) is con- 
sistent with the results from e + e~ and hadroproduction data. 

From the inclusive D° and D* + cross sections a charm production cross section 
in deep inelastic ep scattering for 10 GeV 2 < Q 2 < 100 GeV 2 and 0.01 < y < 0.7 
of a (ep — > ccX) = (17.4 ± 1.6 ± 1.7 ± 1.4) nb has been derived. The data have been 
compared to different NLO calculations including the HI QCD fit to the total F 2 data. 
The charm production cross section is found to be somewhat larger than predicted. 

A first measurement of the charm contribution F 2 ° [x, Q 2 ) to the proton structure 
function for Bjorken x between 8 • 10 -4 and 8 • 10 -3 has been performed. Comparison of 
the present result with the EMC data reveals a steep rise of F 2 C (x, Q 2 ) with decreasing 
x. Agreement is observed between the measured F 2 C and the result of the NLO QCD fit 
of HI to the inclusive F 2 data. Averaged over the kinematic range a ratio I^F 2 C / F 2 ^ = 

0.237 ± 0.021 t° ° tl is obtained, which is one order of magnitude larger than at larger 
x. 
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